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fvq Abstract 

'IIj' Consider a set of discounted optimal stopping problems for a one-parameter family of 

^^ objective functions and a fixed diffusion process, started at a fixed point. A standard problem 

^H in stochastic control/optimal stopping is to solve for the problem value in this setting. 

^SJ In this article we consider an inverse problem; given the set of problem values for a family 

of objective functions, we aim to recover the diffusion. Under a natural assumption on the 

family of objective functions we can characterise existence and uniqueness of a diffusion for 

which the optimal stopping problems have the specified values. The solution of the problem 

relics on techniques from generalised convexity theory. 

Keywords: optimal stopping, generalised convexity, generalised diffusions, inverse American 
3 option problem 

^— I 1 Introduction 
> 

O 

^^ Consider a classical optimal stopping problem in which we are given a discount parameter, an 

^^ objective function and a time-homogeneous diffusion process started at a fixed point, and we 

yyj are asked to maximise the expected discounted payoff. Here the payoff is the objective function 

(^ evaluated at the value of the diffusion at a suitably chosen stopping time. We call this problem 

^^ the forward optimal stopping problem, and the expected payoff under the optimal stopping rule 

. . the (forward) problem value. 






X 



The set-up can be generalised to a one-parameter family of objective functions to give a one- 
parameter family of problem values. In this article we are interested in an associated inverse 
C^ problem. The inverse problem is, given a one-parameter family of objective functions and 

associated optimal values, to recover the underlying diffusion, or family of diffusions, for which 
the family of forward stopping problems yield the given values. 

The approach of this article is to exploit the structure of the optimal control problem and the 
theory of generalised convexity from convex analysis to obtain a duality relation between the 
Laplace transform of the first hitting time and the set of problem values. The Laplace transform 
can then be inverted to give the diffusion process. 
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The generalised convexity approach sets this article apart from previous work on this problem, 
see pi m [6] . All these papers are set in the realm of mathematical finance where the values of 
the stopping problems can be identified with the prices of perpetual American options, and the 
diffusion process is the underlying stock process. In that context, it is a natural question to ask: 
Given a set of perpetual American option prices from the market, parameterised by the strike, 
is it possible to identify a model consistent with all those prices simultaneously? In this article 
we abstract from the finance setting and ask a more general question: When can we identify a 
time-homogeneous diffusion for which the values of a parameterised family of optimal stopping 
problems coincide with a pre-specified function of the parameter. 

Under restrictive smoothness assumptions on the volatility coefficients, Alfonsi and Jourdain [2] 
develop a 'put-call parity' which relates the prices of perpetual American puts (as a function of 
strike) under one model to the prices of perpetual American calls (as a function of the initial 
value of the underlying asset) under another model. This correspondence is extended to other 
payoffs in [1]. The result is then applied to solve the inverse problem described above. In both 
papers the idea is to find a coupled pair of free-boundary problems, the solutions of which can 
be used to give a relationship between the pair of model volatilities. 

In contrast, in Ekstrom and Hobson |6j the idea is to solve the inverse problem by exploiting a 
duality between the put price and the Laplace transform of the first hitting time. This duality 
gives a direct approach to the inverse problem. It is based on a convex duality which requires 
no smoothness on the volatilities or option prices. 

In this article we consider a general inverse problem of how to recover a diffusion which is consis- 
tent with a given set of values for a family of optimal stopping problems. The solution requires 
the use of generalised, or n-convexity (Carlier [1], Villani [H], Rachev and Riischendorf |11]). 
The log- value function is the u-convex dual of the log-eigenfunction of the generator (and vice- 
versa) and the u-subdifferential corresponds to the optimal stopping threshold. These simple 
concepts give a direct and probabilistic approach to the inverse problem which contrasts with 
the involved calculations in [21 [1] in which pdes play a key role. 

A major advantage of the dual approach is that there are no smoothness conditions on the value 
function or on the diffusion. In particular, it is convenient to work with generalised diffusions 
which are specified by the speed measure (which may have atoms, and intervals which have zero 
mass) . 

Acknowledgement: DGH would like to thank Nizar Touzi for suggesting generalised convexity 
as an approach for this problem. 



2 The Forward and the Inverse Problems 

Let Af be a class of diffusion processes, let p be a discount parameter, and let Q = {G{x, 6); 6 £ 
0} be a family of non-negative objective functions, parameterised by a real parameter 6 which 
lies in an interval 0. The forward problem, which is standard in optimal stopping, is for a given 
X G Af, to calculate for each 9 G Q, the problem value 

V{e) = Vx{e) = snpEo[e-'^G{Xr,9)], (2.1) 

T 

where the supremum is taken over finite stopping times r, and Eo denotes the fact that Xq = 0. 
The inverse problem is, given a fixed p and the family Q, to determine whether V = {V{9) : 6 G 



0} could have arisen as a solution to the family of problems (2.1 ) and if so, to characterise those 



elements X £ X which would lead to the value function V. The inverse problem, which is the 
main object of our analysis, is much less standard than the forward problem, but has recently 
been the subject of some studies ([2l[Tl[6]) in the context of perpetual American options. In these 
papers the space of candidate diffusions is Xgtock-, where Xstock is the set of price processes which, 
when discounted, are martingales and G{x, 9) = [0 — x)'^ is the put option payoff (slightly more 
general payoffs are considered in pQ). The aim is to find a stochastic model which is consistent 
with an observed continuum of perpetual put prices. 

In fact it will be convenient in this article to extend the set X to include the set of generalised 
diffusions in the sense of Ito and McKean [8j. These diffusions are generalised in the sense 
that the speed measure may include atoms, or regions with zero or infinite mass. Generalised 



diffusions can be constructed as time changes of Brownian Motion, see Section 5.1 below, and 
also ^, [TU], [13], and for a setup related to the one considered here, [H]. 

We will concentrate on the set of generalised diffusions started and reflected at 0, which are 
local martingales (at least when away from zero). We denote this class Xq. (Alternatively we 
can think of an element X as the modulus of a local martingale Y whose characteristics are 
symmetric about the initial point zero.) The twin reasons for focusing on Xq rather than X , are 
that the optimal stopping problem is guaranteed to become one-sided rather than two-sided, and 
that within Xq there is some hope of finding a unique solution to the inverse problem. The former 
reason is more fundamental (we will comment in Section |6.2| below on other plausible choices 
of subsets of X for which a similar approach is equally fruitful) . For X G Aq , is a refiecting 
boundary and we assume a natural right boundary but we do not exclude the possibility that 
it is absorbing. Away from zero the process is in natural scale and can be characterised by its 
speed measure, and in the case of a classical diffusion by the diffusion coefficient a. In that case 
we may consider X G Aq to be a solution of the SDE (with reflection) 

dXt = a{Xt)dBt + dLt Xq = 0, 

where L is the local time at zero. 

We return to the (forward) optimal stopping problem: For flxed X define ip{x) = ipx{x) = 
Eo[e~''^'']~^, where Hx is the first hitting time of level x. Let 



V{9)= sup [G{x,e)^Q[e-P"^]]= sup 

x:<fi{x)<oo x:ip(x)<oo 



G{x, 



(2.2) 



ip{x) 
Clearly V > V. Indeed, as the following lemma shows, there is equality and for the forward 



problem (2.1), the search over all stopping times can be reduced to a search over first hitting 
times. 

Lemma 2.1. V and V coincide. 

Proof. See Appendix. D 

The first step in our approach will be to take logarithms which converts a multiplicative problem 
into an additive one. Introduce the notation 



log(p{x). 



v{0) -- 


= log(^(^)), 


gix,9) -- 


= log{G{x,9)) 


V(x) = 


= log(Eo[e-^^- 



Then the equivalent log-transformed problem (compare ( |2.2[ )) is 

vie) = snp[g{x,e)-ip{x)], (2.3) 

X 

where the supremum is taken over those x for which ^(x) is finite. To each of these quantities 
we may attach the superscript X if we wish to associate the solution of the forward problem to 



a particular diffusion. For reasons which will become apparent, see Equation (2.5) below, we 
call ifx the eigenfunction (and ipx the log-eigenfunction) associated with X. 

In the case where g{x, 9) = Ox, v and ip are convex duals. More generally the relationship 
between v and tp is that of n-convexity ([1], |14j . |11]). (In Section [3] we give the definition of 
the M-convex dual /" of a function /, and derive those properties that we will need.) For our 
setting, and under mild regularity assumptions on the functions g, see Assumption |3.6| below, 
we will show that there is a duality relation between v and ^p via the log-payoff function g which 
can be exploited to solve both the forward and inverse problems. In particular our main results 



(see Proposition 4.4 and Theorems 5.1 and 5.4 for precise statements) include: 



Forvifard Problem: Given a diffusion X e Xq, let ipxix) = {Eoi^'''^"'])^^ and ipx{x) = 
log{ipx{x)). Set ip^{6) = sup^{g{x,6) — Tp{x)}. Then the solution to the forward problem is 
given by V{6) = exp^ip^^O)), at least for those 6 for which there is an optimal, finite stopping 
rule. We also find that V is locally Lipschitz over the same range of 6. 



Inverse Problem: For v = {v{9) : 9 G Q = [9-, 9^]} to be logarithm of the solution of (2.1) 
for some X G Afo it is sufficient that the g-convex dual (given by v^{x) = supg{g{x, 9) — v{9)}) 
satisfies v^{0) = 0, e^^^^' is convex and increasing, and v^{x) > {g{x,9-) — g{0,9^)} for all 
x>0. 

Note that in stating the result for the inverse problem we have assumed that O contains its 
endpoints, but this is not necessary, and our theory will allow for Q to be open and/or unbounded 
at either end. 

If X is a solution of the inverse problem then we will say that X is consistent with {V{9);9 £ 0}. 
By abuse of notation we will say that ipx (or tpx) is consistent with V (or v = logF) if, when 



solving the optimal stopping problem (2.1) for the diffusion with eigenfunction ipx, we obtain 
the problem values V{9) for each 9 £ Q. 



The main technique in the proofs of these results is to exploit (2.3) to relate the fundamental 
solution ip with V. Then there is a second part of the problem which is to relate ^ to an element 
of X. In the case where we restrict attention to Xq, each increasing convex ip with (p{0) = 1 is 
associated with a unique generalised diffusion X £ Xq. Other choices of subclasses of X may or 



may not have this uniqueness property. See the discussion in Section |5.6 
The following examples give an idea of the scope of the problem: 

Example 2.2. Forward Problem: Suppose G{x,9) = e^^ . Let m> 1 and suppose that X £ Xq 
solves dX = a{X)dW + dL for ^(x)-^ = (x'^i"^-^) + (m - l)x™-2)/(2/j). For such a diffusion 

-I -, TTT. 

c^(x) = exp(^x™),x>0. Thenfor9£Q = (0,oo), V{9) = exp(^0^^^). 

Example 2.3. Forward Problem: Let X be reflecting Brownian Motion on the positive half-line 
with a natural boundary at oo. Then (p{x) = cosh.(x^/2p). Let g(x,9) = 9x so that g-convexity 
is standard convexity, and suppose Q = (0,oo). Then 



v{9) = sup[0x — log(cosh(xy 2/>))]. 



It is easy to ascertain that the supremum is attained at x = x*{9) where 



x*(^) = -^tanh-i(-^) (2.4) 



for e G [0, v^). Hence, for 9 £ (0, ^2^) 






v{9) = tanh I 1 — log I coshtanh 



2p V\/2py V VV2p 



wit/i limits u (0) = and v{^/2f)) = log 2. For ^ > ^/2p we have v{9) = oo. 

Example 2.4. Inverse Problem: Suppose that g{x,9) = 9x and @ = (0,^/2p). Suppose also 
that for9Ge 



Then X is reflecting Brownian Motion. 

Note that X £ Xq is uniquely determined, and its diffusion coefficient is specified on M+. In 
particular, if we expand the domain of definition of to (0, oo) then for consistency we must 
have V{Q) = oo for 9 > ^/2p. 
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Example 2.5. Inverse Problem: Suppose G{x,9) = x^ and V{9) = { — ^ "^ ' : 9 G (1,2)}. 

Then ip{x) = l + x"^ for x > 1 and, at least whilst Xt> 1, X solves the SDE dX = p{\ + X)'^dW . 
In particular, V does not contain enough information to determine a unique consistent diffusion 
in Xq since there is some indeterminacy of the diffusion co-efficient on (0, 1). 

Example 2.6. Inverse Problem: Suppose g{x,9) = —9^/{2{l + x}), Q = [l,oo) and v{9) = 
{—1/2 — log^ : ^ > 1}. Then the g-dual of v is given by v^{x) = log(l + x)/2, x > and is a 
candidate for ip. However e"^^^' = \Jl + x is not convex. There is no diffusion in Xq consistent 
with V. 

Example 2.7. Forward and Inverse Problem: In special cases, the optimal strategy in the for- 
ward problem may be to 'stop at the first hitting time of infinity' or to 'wait forever'. Nonetheless, 
it is possible to solve the forward and inverse problems. 

Let h be an increasing, differentiate function on [0, oo) with h{0) = 1, such that e^ is convex; let 
f be a positive, increasing, differentiable function on [0, oo) such that \un.x^oo f{x) = 1; and let 
w{9) be a non-negative, increasing and differentiable function on Q = [9-, 9+] with w{9-) = 0. 

Suppose that 

g{x,9) = h{x) + f{x)w{9). 

Note that the cross-derivative gxeix,9) = f'{x)w'{9) is non-negative. 

Consider the forward problem. Suppose we are given a diffusion in Xq with log- eigenf unction 
ip = h. Then the log-problem value v is given by 

v{9) = ip3{9) = sup{g{x,9)-ip{x)} = limsup{/(x)u;(^)} = w{9). 

Conversely, suppose we are given the value function V = e'^ on Q. Then 

w%x) = sup{g(x,6l) - w{9)} = sup{/i(x) + {f{x) - l)w{9)} = h{x) 

is the log- eigenf unction of a diffusion X £ Xq which solves the inverse problem. 



A generalised diffusion X £ Xq can be identified by its speed measure m. Let m he a, non- 
negative, non-decreasing and right-continuous function which defines a measure on M^, and let 
m be identically zero on M~. We call x a point of growth of m if m{xi) < m{x2) whenever 
j;i < X < X2 and denote the closed set of points of growth by E. Then m may assign mass to 
or not, but in either case we assume £ E. We also assume that if ^ = sup{x : x € E} then 
^ + m{^+) = oo. If ^ < oo then either ^ is an absorbing endpoint, or X does not reach ^ in 
finite time. 

The diffusion X with speed measure m is defined on [0, ^) and is constructed via a time-change 
of Brownian motion as follows. 

Let F = {J-u)u>o be a filtration supporting a Brownian Motion B started at with a local time 
process {L^;u > 0,z G M}. Define T to be the left-continuous, increasing, additive functional 

Tu= Llm{dz), 

JR 

and define its right-continuous inverse by 

At = inf{n : r„ > t}. 

If we set Xt = B{At) then Xt is a generalised diffusion which is a local martingale away from 
0, and which is absorbed the first time that B hits ^. 

For a given diffusion X £ Xq recall that (p{x) = ipx{x) is defined via tpx{x) = (Eo[e~''''^^])~"^. It 
is well known (see for example [l3l V.50] and [5, pp 147-152]) that (px is the unique increasing, 
convex solution to the differential equation 

j/f = /'/; /(0) = 1, /'(0-) = 0. (2.5) 

2 amax 



Conversely, given an increasing convex function ip with 9^(0) = 1 and <p'{0+) > 0, (2.5) can 
be used to define a measure m which in turn is the speed measure of a generalised diffusion 
X eXo. 

If m{{x}) > then the process X spends a positive amount of time at x. If x £ E is an isolated 
point, then there is a positive holding time at x, conversely, if for each neighbourhood N^ of x, 
m also assigns positive mass to Nx \ {x}, then x is a sticky point. 

li X £ Xq and m has a density, then m{dx) = a{x)~'^dx where a is the diffusion coefficient of 



X and the differential equation (2.5) becomes 



]^a{xff"{x)-pf{x) = G. (2.6) 

In this case, depending on the smoothness of g, v will also inherit smoothness properties. 



Conversely, 'nice' v will be associated with processes solving (2.6) for a smooth a. However 



rather than pursuing issues of regularity, we prefer to work with generalised diffusions. 



3 u-convex Analysis 

In the following we will consider u-convex functions for u = u{y, z) a function of two variables 
y and z. There will be complete symmetry in role between y and z so that although we will 



discuss li-convexity for functions of y, the same ideas apply immediately to n-convexity in the 
variable z. Then, in the sequel we will apply these results for the function g^ and we will apply 
them for (7-convex functions of both x and Q. 

For a more detailed development of u-convexity, see [11] . |14j . [1] and the references therein. 
Proofs of the results below are included in the Appendix. 

Let Dy and D^, be sub-intervals of M. We suppose that u : Dy x D^ 1— )■ M is well defined, though 
possibly infinite valued. 

Definition 3.1. / : Dy — )• M+ is u-convex iff there exists a non-empty S C D^ x M such that 
for all y ^ Dy 

f{y) = sup [u{y,z) + a]. 

{z,a)eS 

Definition 3.2. The u-dual of / is the u-convex function on Dz given by 

nz)= supHy,z)-f{y)]. 
y&Dy 

A fundamental fact from the theory of ti-convexity is the following: 
Lemma 3.3. A function f is u-convex iff (/")" = /. 

The function (/")" (the tt-convexification of /) is the greatest u-convex minorant of / (see the 
Appendix). The condition (/")" = / provides an alternative definition of a n-convex function, 
and is often preferred; checking whether (/")" = / is usually more natural than trying to 
identify the set S. 

Diagrammatically (see Figure 1.), we can think of — (/")(z) = miy[f{y) — u{y, z)] as the vertical 
distance between / and u(., z). Thus /"(z) < when f{y) > u{y, z) for all y £ Dy. 

The following description due to Villani |14j is helpful in visualising what is going on: / is 
li-convex if at every point y we can find a parameter z so that we can caress f from below with 

u{.,z). 

The definition of the u-dual implies a generalised version of the Young inequality (familiar from 
convex analysis, e.g [l2j). 

f{y) + nz)>u{y,z) 

for all {y,z) £ Dy x Dz- Equality holds at pairs {y,z) where the supremum 

sup[u(y,z)-/"(z)] 

z 

is achieved. 

Definition 3.4. The u-subdifferential of f at y is defined by 

5"/(y) = {z€Dz: f{y) + f^{z) = u{y, z)}, 
or equivalently 

5"/(y) = {zeDz: u{y, z) - f{y) > n(y, z) - f{y),yy G Dy}. 

If f/ is a subset of Dy then we define d^f{U) to be the union of n-subdifferentials of / over all 
points in U. 




Figure 1: / is n-subdifferentiable. i9"/(yi) = z*{yi) and 9"(y2) = z*{y2) for 2/2 £ (2/2' ^2)- '^^^ 
distance between m(.,2;) and / is equal to -/'""{z). Note that the u-subdifferential is constant 
over the interval (2/2,2/2) 



Definition 3.5. / is u-subdifferentiable at y if d^f{y) 7^ 0. /is u-subdifferentiable on U if it 
is u-subdifferentiable for all y & U, and / is u-subdifferentiable if it is n-subdifferentiable on 



U 



Dy. 



In what follows it will be assumed that the function u{y, z) is satisfies the following 'regularity 
conditions'. 

Assumption 3.6. (a) u{y^ z) is continuously twice differentiable. 



(b) Uy{y,z) = ^u{y,z) as a function of z, and Uz{y,z) 
strictly increasing. 



;u{y, z) as a function of y, are 



Remark 3.7. We will see below that by assuming |3.6[ a) irregularities in the value function (2.1 ) 
can be identified with extremal behaviour of the diffusion. 



Remark 3.8. Condition |3.6[ b) is known as the single crossing property and as the Spence- 
Mirrlees condition ([!]). If instead we have the 'Reverse Spence-Mirrlees condition': 
(bb) Uy{y, z) as a function of z, and Uz{y, z) as a function of y, are strictly decreasing, 
then there is a parallel theory, see Remark |3. 12 



The following results from n-convex analysis will be fundamental in our application of n-convex 
analysis to finding the solutions of the forward and inverse problems. 



Lemma 3.9. Suppose f is u-subdifferentiable, and u satisfies Assumption 3.6. Then d^f is 

monotone in the following sense: 

Let y,y £ Dy, y > y. Suppose z G d^f{y) and z G d^f{y). Then z > z. 

Definition 3.10. We say that a function is strictly n-convex, when its n-subdifferential is 
strictly monotone. 



Proposition 3.11. Suppose that u satisfies Assumption 3.6. 



Suppose f is a.e differentiable and u-suhdifferentiahle. Then there exists a map z* : Dy — t- D^ 
such that if f is differentiable at y then f{y) = u{y,z*{y)) — f^{z*{y)) and 

f'{y) = uy{y,z*{y)). (3.1) 

Moreover, z* is such that z*{y) is non- decreasing. 

Conversely, suppose that f is a.e differentiable and equal to the integral of its derivative. If 
^3.1\ ) holds for a non- decreasing function z*{y), then f is u-convex and u-subdifferentiable with 
f{y)=u{y,z*{y))-f^{z*{y)). 

Note that the subdifferential d^ f{y) may be an interval in which case z*{y) may be taken to be 



any element in that interval. Under Assumption 3.6, z*{y) is non-decreasing 



We observe that since u{y,z*{y)) = f{y) + f^{z*{y)) we have u{y*{z),z) = f{y*{z)) + f{z) 
and y*{z) E d^f^{z) so that y* may be defined directly as an element of 5"/". If z* is strictly 
increasing then y* is just the inverse of z* . 



Remark 3.12. If n satisfies the 'Reverse Spence-Mirrlees' condition, the conclusion of Lemma 3.9 



is unchanged except that now 'z > z\ Similarly, Proposition 3.11 remains true, except that 
z*{y) and y*{z) are non-increasing. 



Proposition 3.13. Suppose that u satisfies Assumption 3.6 



Suppose f is u-subdifferentiable in a neighbourhood of y. Then f is continuously differentiable 
at y if and only if z* is continuous at y. 



4 Application of w-convex analysis to the Forward Problems 



Now we return to the context of the family of optimal control problems (2.1 ) and the represen- 



tation (2.3). 



Lemma 4.1. Let X £ Xq be a diffusion in natural scale reflected at the origin with a finite or 
infinite right boundary point (,. Then the increasing log-eigenfunction of the generator 

V'x(:E) = -log(E[e-^^1-') 
is locally Lipschitz continuous on (0,^). 

Proof. ipx{x) is increasing, convex and finite and therefore locally Lipschitz on (0,^). 9j(0) = 1, 
and since log is locally Lipschitz on [l,oo), V = log(c/5) is locally Lipschitz on (0,.^). D 



Henceforth we assume that g satisfies Assumption 3.6, so that g is twice differentiable and 
satisfies the Spence-Mirrlees condition. We assume further that G{x,6) is non-decreasing in x. 
Note that this is without loss of generality since it can never be optimal to stop at x' > j; if 
G{x',9) < G{x,9), since to wait until the first hitting time of x' involves greater discounting 
and a lower payoff. 



Consider the forward problem. Suppose the aim is to solve (2.3) for a given X £ Xq with asso- 



ciated log-eigenfunction '([^{x) = ipxi^) = — logEo[e ^^''] for the family of objective functions 



{G{x,9) : 9 G B}. Here is assumed to be an interval with endpoints 0_ and 0+, such that 

Now let 

v{e)= sup [g{x,9)-^{x).] (4.1) 

x:ip(x)<oo 

Then f = -(/'^ is the (^-convex dual of ip. 

By definition d^v{9) = {x : v{9) = g{x, 9) — 'ip{x)} is the (set of) level(s) at which it is optimal to 
stop for the problem parameterised by 9. If d^v{9) is empty then there is no optimal stopping 
strategy in the sense that for any finite stopping rule there is another which involves waiting 
longer and gives a higher problem value. 

Let 9r be the infimum of those values of 6 £ Q such that 

d^v{9) = 0. If v is nowhere /^-subdifferentiable then we set 9ji = 6^. 

Lemma 4.2. The set where v is g-subdijjerentiable forms an interval with endpoints 9- and 
Or. 

Proof. Suppose v is gf-subdifferentiable at 0, and suppose 9 G {9-, 9). We claim that v is 
5-subdifferentiable at 9. 

Fix X G d^v{9). Then v{9) = g{x, 9) — ip{x) and 

g{x,9)-^{x)>g{x,9)-^{x), Vx < e, (4.2) 

and for x = .^ if .^ < oo. We write the remainder of the proof as if we are in the case .^ < oo; 
the case .^ = oo involves replacing x < ^ with x < ^. 

Fix 9 < 9. We want to show 

g{x,9)-i^{x)>g{x,9)-i^ix), VxG(x,e], (4.3) 

for then 

sup{g(a;, 9) - ip{x)} = sup{5'(2;, 9) - ^'(a;)}, 

x<^ x<x 

and since g{x, 9) — %Ij{x) is continuous in x the supremum is attained. 
By assumption, ge{x., t) is increasing in x, and so for x G (x, ^] 



Je 



[ge{x,t) - g0{x,t)]dt < 
Je 

or equivalently, 



g{x,9)-g{x,9)<g{x,9)-g{x,9). (4.4) 



Subtracting (4.4) from (4.2) gives (^4.3|. D 



Lemma 4.3. v is locally Lipschitz on {9^,9ii). 

Proof. On {9 ^,9 ft) v{9) is (7-convex, (7-subdifferentiable and x*{9) is monotone increasing. 

Fix 9', 9" such that 9^ < 9" < 9' < 9r. Choose x' G d9v{9') and x" G d3v{9") and suppose g 
has Lipschitz constant K' (with respect to 9) in a neighbourhood of {x',9'). 

10 



Then v{e') = g{x', 6') - ^{x') and v{e") > g{x', 9") - V'(x') so that 

v{e') - v{e") < g{x', e') - g{x', 9") < K'{9' - 9") 
and a reverse inequahty fohows from considering v{9") = g{x" , 9") — ^(x"). D 



Note that it is not possible under our assumptions to date {g satisfying Assumption 3.6, and g 
monotonic in x) to conclude that v is continuous at 0_, or even that v{9-) exists. Monotonicity 
guarantees that even if 0_ ^ G we can still define x*(6'_) := linig^Q_ x*{9). For example, suppose 
= (0, oo) and for e E (0, 1) let ge{x,9) = g{x,0) + ef{6). Then if v^ij)) is the (^^-convex dual 
of ijj we have Ve{9) = v{9) + ef{9), where v{9) = vo{9). If g and ip are such that lini0^Qv{9) 
exists and is finite, then choosing any bounded / for which limg^o /(^) does not exist gives an 
example for which limg^o^e(^) does not exist. It is even easier to construct modified examples 
such that v{9-) is infinite. 

Denote ^{9,C) = limsup^^^{g{x,9) - iP{x)}. Then for 9r<9 < 9+, ^^9(9) = S(0,^). We have 
shown: 



Proposition 4.4. If g satisfies Assumption 3.6, g is increasing in x and if X is a reflecting 



diffusion in natural scale then the solution to the forward problem is V{9) = exp(V'^(^)). 

Remark 4.5. Suppose now that gx{x, .) is strictly decreasing (the reverse Spence-Mirrlees con- 
dition). The arguments above apply with the obvious modifications. Let 9l be the supremum 



of those values 9^0 such that x*{9) = 0. Then the analogues to Lemmas 4.2 and 4.3 show 



that V is g-subdifferentiable and locally Lipschitz on {9l, 0+) and that for 9^ < 9 < 9i 

V{9) = exp{J:{9,0). 

We close this section with some examples. 



Example 4.6. Recall Example \2.^ hut note that in that example 9 was restricted to take values 
in Q = (1, 2). Suppose Q = [0, oo), g(x, 9) = 9\ogx and 'ijj{x) = log(l + x^). Then 9r = 2 and 
for 9 <9r, x*{9) = (61/(2 - 9))^/'^. Further, for 9 < 2 

v{9) = ^ log{9) + ^ log(2 -9)- log 2, 

and v{9) = oo for 9 > 2. 

Note that v is continuous on [0,9r], but not on 0. 

Example 4.7. Suppose g{x,9) = x9 and = (0,oo). Suppose X is a diffusion on [0, 1), with 

feS^- Thenip{x) = j^ 



1 a natural boundary and diffusion coefficient a{x)^ = ,,^2 ■ Then ip{x) = t^-^ ^'^'^ 



v{9) = sup[6lx + log(l - x^)]. 

X<1 

R is straightforward to calculate that x*{9) = y\ + W^ — 1/9 and then that v{9) : (0, 00) — t- M 
is given by 

v{9) = Vl + 0^-1 -log ( ^J= ). (4.5) 
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5 Application of w-convex analysis to the Inverse Problem 

Given an interval ©CM with endpoints 9- and 0+ and a value function V defined on G we 
now discuss how to determine whether or not there exists a diffusion in Xq that solves the 



inverse problem for V . Theorem 5.1 gives a necessary and sufficient condition for existence. 



This condition is rather indirect, so in Theorem 5.4 we give some sufficient conditions in terms 
of the g-convex dual v^ and associated objects. 

Then, given existence, a supplementary question is whether {V{9) : 6 € 0} contains enough 
information to determine the diffusion uniquely. In Sections |5.3[ |5.4| and |5.5| we consider three 



different phenomena which lead to non-uniqueness. Finally in Section 5.6 we give a simple 
sufficient condition for uniqueness. 

Two key quantities in this section are the lower and upper bound for the range of the g- 
subdifferential of v on 0. Recall that we are assuming that the Spence-Mirrlees condition holds 
so that X* is increasing on Q. Then, if v is somewhere 5f-subdifferentiable we set x- = sup{x G 
d^v{6^)}, or if 9- ^ ©, x_ = \\in.oi0_ x*{9). Similarly, we define 3;+ = inf{x G 9^f(0+)}, or if 
0+ ^ 0, x+ = YmiQ^Q^ x*{9), and xr = lim.Q^Qj^x*{9). If v is nowhere ^-subdifferentiable then 
we set x_ = Xij = x+ = c«. 

5.1 Existence 

In the following we assume that v is g-convex on M-|- x 0, which means that for all G 0, 

v{9) = v3^{9) = sn^{g{x,9) - v^^x)}. 

Trivially this is a necessary condition for the existence of a diffusion such that the solution 
of the optimal stopping problems are given by V . Recall that we are also assuming that g is 



increasing in x and that it satisfies Assumption 3.6 



The following fundamental theorem provides necessary and sufficient conditions for existence of 
a consistent diffusion. 

Theorem 5.1. There exists X ^ Xq such that Vx = V if and only if there exists </> : [0, 00) — t- 
[l,co] such that 0(0) = 1, (j) is increasing and convex and (j) is such that {\og(j))^ = v on Q. 



Proof. If X G Aq then </>x(0) = 1 and 4>x is increasing and convex. Set tpx = log 0x- If Vx = V 
then 

v{9) = vx{0) = sup{5r(x,6') -%Ijx{x)} = ip'^x- 

X 

Conversely, suppose (j) satisfies the conditions of the theorem, and set '0 = log 0. Let £^ = sup{x : 
4>{x) < 00}. Note that if ^ < 00 then 

(log 0)^^) = sup{g{x, 9) - V'(x)} = sup{g{x, 9) - 0(x)} 

x>0 x<^ 

and the maximiser x*{9) satisfies x*{9) < ^. 
For < X < .^ define a measure m via 

, , , 1 0"(x) . r{x) + {^'{x))\ 

m[dx) = ■ — dx = dx. (5.1) 

2p (j) 2/9 
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Let m{dx) = oo for 2; < 0, and, if ^ is finite m{dx) = 00 for x > $,. We interpret (5.1) in a 
distributional sense whenever (j) has a discontinuous derivative. In the language of strings ^ is 
the length of the string with mass distribution m. We assume that ^ > 0. The case ^ = is a 
degenerate case which can be covered by a direct argument. 

Let i? be a Brownian motion started at with local time process LJ and define (r„)„>o via 

JvL Jo 2p i^{Bs) 

Let A be the right-continuous inverse to T. Now set Xt = Bai- Then X is a local martingale 
(whilst away from zero) such that d{X)t/dt = dAt/dt = {dm / dx\x=Xt)^^ ■ When m{dx) = 
a{x)~'^dx, we have d{X)t = a{Xt)'^dt. 

We want to conclude that E[e~''^^] = ex.p{—ip{x)). Now, ^px{x) = {'Kle^^^'-'])^^ is the unique 
increasing solution to 

2 dmdx 
with the boundary conditions /'(O— ) = and /(O) = 1. Equivalently, for all x,y €z {0,(,) with 
X < y, tpx solves 

f\y-)-f\x-)= f 2pf{z)m{dz). 

By the definition of m above it is easily verified that exp('i/'(x)) is a solution to this equation. 
Hence (j) = (px and our candidate process solves the inverse problem. D 

Remark 5.2. Since v is ^r-convex a natural candidate for cj) is e^*^^^, at least if v^{G) = and e"^ 
is convex. Then (j) is the eigenfunction ipx of a diffusion X ^ Xq. 

Our next example is one where 0(a;) = e" '^^ is convex but not twice differentiable, and in 
consequence the consistent diffusion has a sticky point. This illustrates the need to work with 
generalised diffusions. For related examples in a different context see Ekstrom and Hobson [6j. 

Example 5.3. Let = M_|_ and let the objective function he g{x,9) = exp{9x). Suppose 

exp(ie2) < < 2, 

V{e) = { exp(0 - 1) 2 < < 3, 

exp(^e3/2) 3 < ^_ 



Writing ip = e" we calculate 



^i^l - ^ ™^^^3 



exp(x^) < X < 1, 
exp(x'^) 1 < X. 



Note that ip is increasing and convex, and ^{0) = 1. Then ip' jumps at 1 and since 

^{l) = ^'(1+) - ip'il-) = 2p^(l)m({l}) 

we conclude that m({l}) = ^. Then Tu includes a multiple of the local time at 1 and the 
diffusion X is sticky there. 



Theorem 5.1 converts a question about existence of a consistent diffusion into a question about 



existence of a log-eigenfunction with particular properties including {\og4>y = v. We would 
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like to have conditions which apply more directly to the value function V{-). The conditions we 
derive depend on the value of x-. 



As stated in Remark 5.2, a natural candidate for (j) is e^^^^' . As we prove below, if x_ =0 this 
candidate leads to a consistent diffusion provided v^{0) = and e"^^^' is convex and strictly 
increasing. If a;_ > then the sufficient conditions are slightly different, and e" need not be 
globally convex. 

Theorem 5.4. Assume v is g-convex. Each of the following is a sufficient condition for there 
to exist a consistent diffusion: 

1. X- = 0, v^{0) = and e^"'^' is convex and increasing on [0,a;-|_). 

2. < X- < oo, v^{x-) > 0, e""^^^' is convex and increasing on [x_,x+), and on [0,X-), 
v^{x) < f{x) = \og{F{x)) where 

^, , ^ exp(t;5(x_)) — 1 
F(x) = 1 + x^-^ — ^ — 

X- 

is the straight line connecting the points (0, 1) and {x-,e^^^^~'). 

3. X- = oo and there exists a convex, increasing function F with log(F(0)) = such that 
f{x) > v^{x) for all x > and 

lim{f{x)-v3{x) = 0}, 
where f = logF. 



Proof. We treat each of the conditions in turn. If x- =0 then Theorem 5.1 applies directly 
on taking (j){x) = e^^^^', with (j){x) = oo for x > x^ (we use the fact that v is (^-convex and so 
v39 = v). 

Suppose < a;_ < cxd. The condition e^^^^' < F{x) on [0,x„) implies F'{x) = (e""'^-) — 
^)/ X- < {e"" '^- "•*)'. Although the left-derivative v ^{x— )' need not equal the right-derivative 



v^{x+) the arguments in the proof of Proposition 3.11 show that v^{x—) < v^{x+y. This 
implies that the function 

, , . f F(x) X < X- 

0F[X) = < / ru W 

[ exp[v''[x)) X- < X < x+ 

is convex at x^ and hence convex and increasing on [0,2;+). 

Setting (j)p{x+) = lima;f-a;_^ (Pf{x) and (/li? = oo for x > x+ we have a candidate for the function 
in Theorem 15. 1[ 



It remains to show that (log (J)f)^ = u on 0. We now check that (pp is consistent with V on 0, 
which follows if the g-convex dual of ^ = log{(pF) is equal to v on Q. 

Since ip > v^ we have "0^ ^ v. We aim to prove the reverse inequality. By definition, we have 
for 6) G G 



i;3(e) = ( sup {g{x, 9) - f{x)} V sup {g{x, 6) - v^x)} . (5.2) 

\X<X- J \x-<X<Xj^ I 
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Now fix X G [0,x_). For 9 < Oji we have by the definition of the (7-subdifferential 

g{x*ie),e)-v<^{x*i9))>gix,e)-v^x). 

Hence v{6) = sup^>Q{g{x,e) - v9{x)} = sup^>^_{g{x,e) - vS{x)} < ^3{Q). 
Similarly, \i 9 > Or we have for all x' G [0, x_), 

limsupg(x, 9) — v^{x) > g[x' , 9) — v^{x'). 

x—^oo 

and v{9) = limsup^{(7(x, 6*) — v^{x)} = sup^>^_{(7(x, 6*) — v^{x)} < ip^{9). 
Finally, suppose 2;_ = 00. By the definition of /^ and the condition f > v^ we get 

f^i9) = sup{g{x,9)-f{x)} 

x>0 

< sup{g{x,9)-v>^{x)} 

x>0 

= v{9). 

On the other hand 

v{9) = liiaisup{g{x,0) — f{x) + f{x) — v^{x)} 

x—^00 

< limsup{5(x,0)-/(x)}+ lim{f{x)-v3{x)} < f9{9). 

Hence v{9) = f9{0) on 9. D 

Remark 5.5. Case 1 of the Theorem gives the sufficient condition mentioned in the paragraph 
headed Inverse Problem in Section [2l If 0_ G then j;_ = if and only if for all x > 0, 
g{x, 0-) — v^{x) < g{0, 0_), where we use the fact that, by supposition, v^{0) = 0. 

5.2 Non-Uniqueness 

Given existence of a diffusion X which is consistent with the values V{9), the aim of the next 
few sections is to determine whether such a diffusion is unique. 

Fundamentally, there are two natural ways in which uniqueness may fail. Firstly, the domain 
Q may be too small (in extreme cases Q might contain a single element). Roughly speaking the 
g'-convex duality is only sufficient to determine v^ (and hence the candidate (p) over (x_,x+) 
and there can be many different convex extensions of 4> to the real line, for each of which 
■0^ = V. Secondly, even when x- = and x+ = 00, if x*{9) is discontinuous then there can be 
circumstances in which there are a multitude of convex functions (j) with (logc/))^ = v. In that 
case, if there are no for which it is optimal to stop in an interval I, then it is only possible to 
infer a limited amount about the speed measure of the diffusion over that interval. 

In the following lemma we do not assume that ip is ^(-convex. 

Lemma 5.6. Suppose v is g-convex and ipd = v on 0. Let A{9) = {x : g{x, 9) —ip{x) = 11^3(9)}. 
Then, for each 9, A{9) C 5^0^(0) = d3v{9), and for x e A{9), ^{x) = ^39^x) = v3{x). Further, 
for 9 G {9-,9r) we have A{9) / 0. 
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Proof. Note that if ip is any function, with ips = y then ^ > ip^^ = v^. 

If X G A{9) then 

^s(6l) = g{x, 0) - tpix) < g{x, 9) - vHx) < v{e). 

Hence there is equahty throughout, so x G d^v{6) and iplx) = v^(x) = ip^^{x). 

For the final part, suppose < Or and fix G {9,9^. From the Spence-Mirrlees condition, if 

X > X := x*{9), 

g{x, 9) - g{x, 9) < g{x, 9) - g{S:, 9), 

and hence 

{g{x, 9) - ^{x)} - {<7(x, 9) - V'(x)} < {g{x, 9) - V'(x)} - {g{x, 9) - ^{x)} < 0. 
In particular, for x > x, ^(x, 9) — ipix) < g{x, 9) — ip{x) and 

supg{x,9) — ^p{x) = sup g{x.,9) — ip{x). 

x>0 0<x<x 

This last supremum is attained so that A{9) is non-empty. 

D 

5.3 Left extensions 

In the case where x_ > and there exists a diffusion consistent with V then it is generally 
possible to construct many diffusions consistent with V . Typically V contains insufficient in- 
formation to characterise the behaviour of the diffusion near zero. 



Suppose that < x_ < oo. Recall the definition of the straight line F from Theorem 5.4 



Lemma 5.7. Suppose that < x_ < oo and that there exists X & Xq consistent with V . 

Suppose that 9^ > 9- and that v^ and is continuous and differentiable to the right at x_. 
Suppose further that x*{9) > x_ for each 9 > 9^. 

Then, unless either v^{x) = /(x) for some x G [0,x_) or (f^)'(x_ +) = /'(x_), there are many 
diffusions consistent with V. 

Proof. Let <j) be the log-eigenfunction of a diffusion X G Xq which is consistent with V 



If 0_ G then f^(x_) = i/;(x_) by Lemma 5.6 Otherwise the same conclusion holds on taking 



limits, since the convexity of (p necessitates continuity of V'- 

Moreover, taking a sequence 0„ | 0_, and using x(0„) > x*(6'.„— ) > x_ we have 

V>(x(g„))-V.(x„) V3{x{9n))-V^{x.) 

ip (x_ +) = hm ——- = hm ——- = (w^) (x_ +) 

Bnie^ x[9n) - x_ 6i,4e- x[9n) - x_ 

In particular, the conditions on v^ translate directly into conditions about <p. 

Since the straight line F is the largest convex function with F{0) = 1 and F{x-) = e'"^^^-' we 
must have (p < F. 

Then if (p{x) = F[x) for some x G (0, x_) or (p'{x- +) = F'(x_), then convexity of (p guarantees 
(p = F on [0,x_]. 
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Otherwise there is a family of convex, increasing with (f)(0) = 1 and such that v^{x) < 
log(/>(x) < F{x) for X < x^ and (f){x) = (p{x) for x > X-. 



For such a (j), then by the arguments of Case 2 of Theorem 5.1 we have (logc/ii?)^ = v and then 
v^ < log(f) < (/"F imphes v > (log(/>)^ > {log(pFy = v. 

Hence each of is the eigenfunction of a diffusion which is consistent with V. D 



Example 5.8. Recall Example 2.5. in which we have x_ = 1, ^^'(l) = 2 and ip{'i) = 2. We can 
extend ip to x & [0, 1) by (for example) drawing the straight line between (0, 1) and (1, 2) (so that 
for X < 1, ip{x) = 1 + x). With this choice the resulting extended function will be convex, thus 
defining a consistent diffusion on M"*". Note that any convex extension of ip (i.e. any function ip 
such that 0(0) = 1 and c^'(0— ) = 0, (p{x) = p{x) for x > 1) solves the inverse problem, (since 
necessarily ip{x) > 2x = e^^'^^' on (0, 1)). The most natural choice is, perhaps, p{x) = 1 + x^ 
for X £ (0,1). 

Our next lemma covers the degenerate case where there is no optimal stopping rule, and for 



all 9 it is never optimal to stop. Nevertheless, as Example 5.10 below shows, the theory of 
u-convexity as developed in the article still applies. 

Lemma 5.9. Suppose x_ = oo, and that there exists a convex increasing function F with 
F{0) = 1 and such that logF{x) > v^{x) and lima;_i.oo{log-F(x) — v^{x)} = 0. 

Suppose that lim2:_j.oo e^^^^V^ exists in (0,oo] and write k = lim^-^-oo e^^^^V^- If n < oo then 
X is the unique diffusion consistent with V if and only if e"^^-^ -* = 1 + kx' for some x' > 
or limsup^-|-(^(l + kx) — e^^*-^^ = 0. If k = oo then there exist uncountably many diffusions 
consistent with V . 



Proof. The first case follows similar reasoning as Lemma 5.7 above. Note that x i— )• 1 + kx is 
the largest convex function F on [0, oo) such that F{0) = 1 and lima;_i.oo —^ = i^- 

If e""^^^ •* = 1 + Kx' for any x' > 0, or if limsup2,+oo(l + kx) — e""^^^' = then there does not exist 
any convex function lying between 1 + kx and e"^^^> on [0, oo). In particular (j){x) = 1 + kx is 
the unique eigenfunction consistent with V. 

Conversely, if e"" lies strictly below the straight line 1 + kx, and if limsup2,^oQ(l + Kx) — e^^^^-' > 
then it is easy to verify that we can find other increasing convex functions with initial value I, 
satisfying the same limit condition and lying between e" and the line. 

In the second case define Fa{x) = F{x) + ax for a > 0. Then since lima;_j.oo e"^^^' /x = oo we 
have 

Hence Fa is the eigenfunction of another diffusion which is consistent with V. We conclude 
that there exist uncountably many consistent diffusions. D 

Example 5.10. Suppose g{x,9) = x^ + ^tanhx and v{9) = 9 on Q = M+. For this example 
we have that v is nowhere g-subdifferentiable and x_ = oo. Then v^{x) = x^ and each of 
p[x) = e^', 

_, , f l + (e-l)x 0< x< 1 

^(") = 1 e^^ 1 < X, 

and ipa{x) = p{x) + ax for any a G M_|_ is an eigenfunction consistent with V . 
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5.4 Right extensions 

The case of x+ < oo is very similar to the case x- > 0, and typicahy if there exists one diffusion 
X £ Xq which is consistent with V, then there exist many such diffusions. Given X consistent 
with V, the idea is to produce modifications of the eigenfunction ipx which agree with ipx on 
[0,x+], but which are different on (x+,oo). 

Lemma 5.11. Suppose x+ < cxd. Suppose there exists a diffusion X £ Xq such that Vx agrees 
with V on Q. Ifv^{xj^) + [v^)'{x^ +) < oo then there are infinitely many diffusions in Xq which 
are consistent with V . 

Proof. It is sufficient to prove that given convex increasing defined on [0,x+) with </)(0) = 1 
and (log(/>)^ = w on 0, then there are many increasing, convex (j) with defined on [0,cxd) with 
(/>(0) = 1 for which (logc/))^ = v. 



The proof is similar to that of Lemma |5.7[ D 

Example 5.12. Let G{x,e) = ex/{e + x), and 6 = (0,oo). 

Consider the forward problem when X is a reflecting Brownian motion, so that the eigenfunction 
is given by (p{x) = cosh(x-^/2p). Suppose p = 1/2. 

Then {g{x,9) — log(cosh2;)} attains its maximum at the solution x = x*{6) to 

x^tanhx 

e=- . (5.3) 

1 — X tanh X 

It follows that X- = but x+ = liuig-^oo x*{(^) = ^ where A is the positive root of C{X) = and 
£(A) = 1- AtanhA. 

Now consider an inverse problem. Let G and be as above, and suppose p = 1/2. Let x*{6) be 



the solution to (5.3) and let v{6) = g{x*{9),6) — log(coshx*(0)). Then the diffusion with speed 
measure m{dx) = dx (reflecting Brownian motion), is an element of Xq which is consistent 
with {V{6) : 9 £ (0,oo)}. However, this solution of the inverse problem is not unique, and any 
convex function ip with ^{x) = coshx for x < X is the eigenfunction of a consistent diffusion. 
To see this note that for x > x+, v^{x) = li'nig-^^{g{x,9) — v{9)} = log(xcosh(x+)/2;+) so that 
any convex ip with (p{x) = coshx for x < X-^ satisfies ^ > e'" . 

Remark 5.13. If x+ +f^(x+) + (f^)'(a;_|_ +) < oo then one admissible choice is to take = oo. 
This was the implicit choice in the proof of Theorem |5.1[ 



Example 5.14. The following example is 'dual' to Example \27i 

Suppose p = 1/2, g{x, 9) = 9x, Q = (0, oo) and v{9) = log(cosh^). Then v^{x) = 2;tanh~"'^(x) + 
2 log(l — x'^), for X < 1. For x > 1 we have that v^ is infinite. Since v is convex, and g-duality 
is convex duality, we conclude that v is g -convex. Moreover, v^ is convex. Setting ip = v^ we 
have that V'(O) = 0, (f = e"^ is convex and tp^ = v^^ = v. Hence ip is associated with a diffusion 
consistent with V, and this diffusion has an absorbing boundary at ^ = 1. 

For this example we have x+ = 1 and v^{x-^) = log 2, but the left- derivative of v^ is infinite 
at x+ and v^ is infinite to the right of x+ . Thus there is a unique diffusion in Xq which is 
consistent with V . 
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5.5 Non-Uniqueness on [x , x+) 

Even if [x_,x+) is the positive real line, then if x*{9) fails to be continuous it is possible that 
there are multiple diffusions consistent with V. 

Lemma 5.15. Suppose there exists a dijfusion X £ Xq which is consistent with {V{9) : 6 G ©}. 

Suppose the g-subdifferential of v is multivalued, or more generally that x*{9) is not continuous 
on Q. Then there exists an interval I C (x_,x+) where the g-suhdifferential of tp = v^ is 
constant, so that 9*{x) = 0, \/x £ I. If G{x,9) = e^^^'' is strictly convex in x on some 
subinterval of Iq of I then the diffusion X is not the unique element of Xq which is consistent 
with V . 

Proof. First note that if x*{9), is continuous then 9* = x*~ is nowhere constant and hence 



strictly monotone and thus ip=v^ is strictly ^r-convex (recall 3.10). 



Suppose G{x) := G{x, 9) is strictly convex on Iq C /. Then we can choose G such that 

• G = G on /^, 

• G is linear on Iq, 

• G is continuous. 
ThenG(a;) > G{x,9). 

By definition we have 

^{x)=g{x,9)-^3{0^ x^l_ 

Then ^px{x) = G{x,9)/V{9) on /, see Figure 2. 

Let <p be given by 

f "Px^x) onig , 

"fix) = < G{x) J 

Then (p is convex and (p > (p, so that they are associated with different elements of Xq. Let 

ip = \\i(p. 

It remains to show that v := ip9 = ifjd = y , It is immediate from ip > ip that ip^ < ip^ . For the 
converse, observe that 



v{9) = {sup{g{x,9)-'iPix)}]\/ {sup{g{x,9)-^Pix)} 
\xeio / \xeiS J 

= sup{g{x,9) -Tpix)} 



sup{g{x,9) -ipix)} 
xeis 



< v{9). 



D 



Example 5.16. Suppose G{x,9) = e^^ , and Q = (0, oo). Suppose that X is such that ip is 
given by 

/ 9 

X <2 



^p{x) 



4 

(x - 1) 2 < x < 3 

-g- + 2 3 < X 
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Figure 2: The dashed hne represents a function e^^^'^'^'^ for some ^ G G and some c. ^p is given 
by the solid hne on /q and the dashed hne on Iq. Then ip = \og(f) has a (7-section over Iq and G 
is convex there. (^, given by the solid hne is another eigenfunction consistent with V . 



It follows that 



v{e) 



< 1 



3^2 



1 < 



Then d^v{l) is multivalued, and there are a family of diffusions X £ Xq which give the same 
value functions as X . 



In particular we can take 

V'(x) 



^ x<2 

log((e2-e)x + 3e-2e2) 2 < j; < 3 

-p- ~r 7T o ^ X 



Then ip^ = v{6) and ip is a log- eigenfunction. 



5.6 Uniqueness of diffusions consistent with V 

Proposition 5.17. Suppose V is such that x*{6) is continuous on Q, with range the positive 
real line. 

Then there exists at most one diffusion X £ Xq consistent with V . 



Proof. The idea is to show that v^{x) is the only function with g-convex dual v. Suppose ip is 
such that ip3 = V on Q. For each x there is a with x*{6) = x, and moreover d^v{6) = {x}. 
Then by Lemma 5.6, A{0) = {x} and ^[x) = v^(x). D 



Recall that we define Or = supg d^v{9) ^ and if Or > 6^ set xr = lime-|-e^ x*{6). 
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Theorem 5.18. Suppose V is such that v is continuously differentiable on {9^, Or) and that 
X-_ = and xj:? = oo. 

Then there exists at most one diffusion X £ Xq consistent with V. 

Proof. The condition on the range of x*{9) translates into the conditions on x_ and xr, so it is 
sufficient to show that x*{6) is continuous at 6* E (6-, Or) if and only if v is differentiable there. 
This follows from Proposition |3. 13 D 



Corollary 5.19. If the conditions of the Theorem hold but either e^^^^' is not convex orv^{0) ^ 
0, then there is no diffusion X £ Xq which is consistent with {V{9),6 G ©}. 



Example 5.20. Recall Example 2.6. For this example we have x*{9) = 9—1, which on 
Q = (l,oo) is continuous and strictly increasing. Then e""^^^' = \/l + x and by the above 
corollary there is no diffusion consistent with v. 

Remark 5.21. More general but less succinct sufficient conditions for uniqueness can be deduced 



from Lemma 5.7 or Lemma 5.11 For example, if < x_ < x+ = oo, but (t'^)'(x-) 
(1 — e~^^^^~')/ x^ then there is at most one X G Xq which is consistent with V. 



6 Further examples and remarks 

6.1 Birth-Death processes 

We now return to Xq and consider the case when E is made up of isolated points only; whence 
X is a birth-death process on points Xn & E indexed by n G No, with associated exponential 
holding times A^. We assume xq = 0, x„ is increasing, and write x^o = lim„ Xn- 

For a birth-death process the transition probabilities are given by 

¥n,n+l{t) = PnXnt + o{t), 
Fn,n-lit) = qnXnt + o{t), 

where of course qn = I — Pn, with po = 1. By our assumption that, away from zero, (Xt)t>o is 
a martingale, we must have Pn = ^"+^~^" , Then we can write x„ = Xn-i + IttT-i ' ■ Let 

X„ + l-X„-l lli = P' 

1 popi...pn-l 
m{Xn) = T . 

An qiq2---qn 



Then it is easy to verify, but see f^, that (2.5) can be expressed in terms of a second-order 
difference operator 



1 



m{xn) 



fjXn+l) - fjXn) _ fjXn) - /(^n-l) 
Xn+1 Xn Xn Xn—1 



pfiXn) = 0, (6.1) 



with boundary conditions /(O) = 1 and /'(O— ) = 0. 



Let M{x) = Y^^ ^^m{x). In the language of strings, the pair (M, [0,Xoo)) is known as the 
Stieltjes String. If Xqo + M (xqo) < oo the string is called regular and Xqo is a regular boundary 
point, while otherwise the string is called singular, in which case we assume that Xqo is natural 
(see Kac [Hj). 
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In this section we consider the call option payoff, G{x, 9) = {x — 6)^ defined for ^ € B = [Oq, oo). 
This objective function is straight-forward to analyse since the (7-duality corresponds to straight 
lines in the original coordinates. It follows that for the forward problem V is decreasing and 
convex in 9. V is easily seen to be piecewise linear. 

Our focus is on the inverse problem. Note that the solution of this problem involves finding 
the space E and the jump rates A„. Suppose that V is decreasing, convex and piecewise linear. 
Let (^n)neNo be a sequence of increasing real valued parameters with ^o < and 0n increasing 
to infinity, and suppose that V has negative slope Si on each interval (^j,0j+i). Then Si is 
increasing in i and 

V{0) = y(0„) + {e- en)Sn for a„ < ^ < dn+l. (6.2) 

We assume that sq = wg-r < 0. 



Since V is convex, v is log((x — 9)'^) convex. Let '^{x) = exp(f^(x)). By Proposition 3.11, for 

9 G [9n,9n+l) 



-1 



ge{x*{9),9) 



:,*^e)-e '"'' ' " ' v{9n) + {9-en)sn 

so that Xn '■= x*{9n) = 9n — V{9n)/sn- Note that x*{0) is constant on [9n,9n+i)- We find that 

fovee[9n,9n+l) 

4^{x*{9)) = log(^„ - 9 - V{9n)/Sn) - v{9), 
and hence ip{x*{9)) = ^. Then, for x G [x*{9n),x*{9n+i)), 

^P{x) = log{x-9n)-v{9n). (6.3) 

We proceed by determining the Q-matrix for the process on [x*{0) = 0,^)- For each n, let pn 
denote the probability of jumping to state Xn+i and g„ the probability of jumping to Xn-i- Then 
Pn and Qn are determined by the martingale property (and po = 1). Further A^ is determined 



either through (6.1) or from a standard recurrence relation for first hitting times of birth-death 
processes: 

Example 6.1. Suppose that 9n = n + 2~^ — 2 so that 9q = —1, and V{9n) = 2"". It follows that 
Sn = — (2"^"*^ — 1)~^. We find Xn = n (this example has been crafted to ensure that the birth-death 
process has the integers as state space, and this is not a general result). Also (p{n) = 2"+^ — 1 
(if is piecewise linear with kinks at the integers) and the holding time at Xn is exponential with 
rate A„ = 4p(l - 2" ("+!)). 

6.2 Subsets of X and uniqueness 

So far in this article we have concentrated on the class Xq. However, the methods and ideas 
translate to other classes of diffusions. 

Let Xm,s denote the set of all diffusions reflected at 0. Here m denotes the speed measure, and 



s the scale function. With the boundary conditions as in (2.5), (p{x) = fx{x) is the increasing, 
but not necessarily convex solution to 

2 dmds 
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In the smooth case, when m has a density m{dx) = v{x)dx and s" is continuous, (6.4) is 
equivalent to 

]^a\x)r{x) + /i(x)/'(x) = pf{x), (6.5) 

where 

u[x) = a-2(x)e^(^), s'{x) = e-*^(^'), M{x) = T 2a-'^ {z) ^i{z)dz , 

Jo- 
see [3]. 

Now suppose V = {V{9) : ^ G 0} is given such that f"(0) = 0, (f")'(0) = and v^ is increasing. 



then we will be able to find several pairs {cr,fi) such that exp(ij'*) solves (6.5) so that there is a 



family of diffusions rather than a unique diffusion in Xm,s consistent with v. 

It is only by considering subsets of Xm,s, such as taking s{x) = as in the majority of this 
article, or perhaps by setting the diffusion co-efficient equal to unity, that we can hope to find 
a unique solution to the inverse problem. 



Example 6.2. Consider Example 2.6 where we found ip{x) = \/l + x. Let Xi^s be the set of 
diffusions with unit variance and scale function s (which are reflected at 0). Then there exists 
a unique diffusion in Xi^s consistent with V . The drift is given by 

^(^,^1/4 + 2,(1+.)= 



1 + x 

7 Applications to Finance 

7.1 Applications to Finance 

Let ^stock be the set of diffusions with the representation 

dXt = {p- 5)Xtdt + r]{Xt)XtdWt. 

In finance this SDE is often used to model a stock price process, with the interpretation that p 
is the interest rate, 6 is the proportional dividend, and r] is the level dependent volatility. Let 
xq denote the starting level of the diffusion and suppose that is an absorbing barrier. 

Our aspiration is to recover the underlying model, assumed to be an element of Xgtock, given a 
set of perpetual American option prices, parameterised by 9. The canonical example is when 



6 is the strike, and G{x,9) = {9 — xy, and then, as discussed in Section 6.2, the fundamental 
ideas pass over from Xq to X stock- We suppose p and 8 are given and aim to recover the volatility 
r/. 

Let ip be the convex and decreasing solution to the differential equation 

\ri{xfx^f.. + {p-5)xf^-pf = Q. (7.1) 

(The fact that we now work with decreasing 99 does not invalidate the method, though it is now 
appropriate to use payoffs G which are monotonic decreasing in x.) Then r/ is determined by 
the Black-Scholes equation 

, = , Mx)-(p-W(^) _ (72) 

x'^ip"[x) 
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Let G = G{x, 6) be a family of payoff functions satisfying assumption 3.6, Under the additional 

assumption that G is decreasing in x (for example, the put payoff) Lemma 2.1 shows that the 

optimal stopping problem reduces to searching over first hitting times of levels x < xq. Suppose 

that {V{9); G 0} is used to determine a smooth, convex ip = exp('i/;) on [0, ^) via the ^'-convex 

duality 

'4^{x) = v^{x) = sup[g{x, 6) — v{6)\. 
See 



Then the inverse problem is solved by the diffusion with volatility given by the solution of (7.2) 
above. Similarly, given a diffusion X G Xstock such that tp = log(c^) is ^r-convex on [0,.^), then 
the value function for the optimal stopping problem is given exactly as in Proposition 4.4. See 
Ekstrom and Hobson [6j for more details. 

Remark 7.1. The irregular case of a generalised diffusion requires the introduction of a scale 
function, see Ekstrom and Hobson [6] . That article restricts itself to the case of the put/call 
payoffs for diffusions Xstock^ using arguments from regular convex analysis to develop the duality 
relation. However, the construction of the scale function is independent of the form of the payoff 
function and is wholly transferable to the setting of (7-convexity used in this paper. 

7.2 The Put-Call Duality; ^f-dual Processes 

In [2], Alfonsi and Jourdain successfully obtain what they term an American Put-Call duality 
(see below) and in doing so, they solve forward and inverse problems for diffusions X G Xgtock- 
They consider objective functions corresponding to calls and puts: G{x,9) = (x — 9)~^ and 
G{x,6) = {9 — x)'^. In [I\ the procedure is generalised slightly to payoff functions sharing 
'global properties' with the put/call payoff. In our notation the assumptions they make are the 
following: [Jj (1.4)] 

Let G : M"*" x M"*" — )■ M"*" be a continuous function such that on the space ^ = {{x,9) : G{x, 9) > 
0} / 0, G is C^ and further for ah x, 6* G $ 

G,(x,0)<O, Ge{x,9)>{), C(x,0)<O Gee(x,^)<0. 

Subsequently, they assume [U (3.4)] 

GG^e > G^Ge on $. (7.3) 



Condition (7.3) is precisely the Spence-Mirrlees condition applied to g = log(G). Note that 
unlike Alfonsi and Jourdain [1] we make no concavity assumptions on G. We also treat the 
case of the reverse Spence-Mirrlees condition, and we consider classes of diffusions other than 
'^stock- Further, as in Ekstrom and Hobson [6] we allow for generalised diffusions. This is 
important if we are to construct solutions to the inverse problem when the value functions are 
not smooth. Moreover, even when v is smooth, if it contains a g-section, the diffusion which is 
consistent with v exists only in the generalised sense. When v contains a g-section we are able to 
address the question of uniqueness. Uniqueness is automatic under the additional monotonicity 
assumptions of [T]. 

In addition to the solution of the inverse problem, a further aim of Alfonsi and Jourdain [2] is 
to construct dual pairs of models such that call prices (thought of as a function of strike) under 
one model become put prices (thought of as a function of the value of the underlying stock) in 
the other. This extends a standard approach in optimal control/mathematical finance where 
pricing problems are solved via a Hamilton-Jacobi-Bellman equation which derives the value of 
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an option for all values of the underlying stock simultaneously, even though only one of those 
values is relevant. Conversely, on a given date, options with several strikes may be traded. 

The calculations given in [2] are most impressive, especially given the complicated and implicit 
nature of the results (see for example, p3 Theorem 3.2]; typically the diffusion coefficients are 
only specified up to the solution of one or more differential equations). In contrast, our results on 
the inverse problem are often fully explicit. We find progress easier since the g-convex method 
does not require the solution of an intermediate differential equation. Furthermore, the main 
duality (pQ Theorem 3.1]) is nothing but the observation that the (7-subdifferentials (respectively 
X* and 6*) are inverses of each other, which is immediate in the generalised convexity framework 
developed in this article. Coupled with the identity, gx{x*{9),9) = 7p'{x*{6)) the main duality 
result of [2] follows after algebraic manipulations, at least in the smooth case where x* and 9* 
are strictly increasing and differentiable, which is the setting of [5] and [T]. 



A Proofs 

A.l The Forward and Inverse Problems 



Proof of Lemma 



2.1 



Clearly V >V, since the supremum over first hitting times must be less 



than or equal to the supremum over all stopping times. 

Conversely, by ( 
local martingale 



Conversely, by (2.2), (p{x) > y'' . Moreover, (2.5) implies that e ^^ip{Xt) is a non-negative 



hus for stopping times r we have 



and hence a supermartingale. Tl 

1 > Eo[e-''XXO] > Eo[e-P^GiXr,9)/V{9)] 
and hence V{9) > sup^ Eo[e-P^G{xr,9)]. 

A. 2 u-convexity 

The methods and many of the results we use in this section are to be found in |llj and [1|. 
Lemma A.l. For every function f : D^ — )• M, (/")" is the largest u-convex niinorant of f. 



Proof. Note that /* = —00 is always a n-convex minorant and by Definition 3.1 a function 
defined as the supremum of a collection of n-convex functions is n-convex. 

Hence we can define /™, the largest n-convex minorant of / by 

r(^)=sup{C(:E);C</,CG^(^.)}, 

C 

where U{Dz) is the set of n-convex functions on Dz- Then by the Fenchel Inequality, / > (/")" 
and since (/")" is n-convex, it follows that /"* > (/")". 

Since /™ is n-convex there exists S' C D^ x M such that 

/""(y) = sup [u{y,z) + a] 

For fixed {z, a) G S", l{y) = u{y, z) + a \s n-convex and it is easy to check that {l^)^{y) = l{y)- 

We note that if / < /, then /" > /". Since u{y,z) + a < f"^{y) < f{y), we therefore have 
u{y, z) + a < (/")"(y). Taking the supremum over pairs {z, a) E S' , we get /"^ < (f^)^- □ 
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Proof of Lemma 3.3 If / is u-convex then / is its greatest tt-convex minorant, so / = (/")'^. 



On the other hand if (/")" = / then since (/")" is u-convex, so is /. 



Proof of Lemma 3.9 Let y < y. For any z £ d^ f{y) and any z £ d'^f{y), the following pair 



of inequalities follow from n-convexity (see especially the second part of Definition 3.4): 



f{y)-f{y) > u{y,z)-u{y,z), 
f{y)-f{y) > u{y,z)-u{y,z). 

Adding and re- arranging we get, 

u{y, z) - u{y, z) + u{y, z) - u{y, z) < 0, 

and hence, 

rv 

/ [uy{v,z) — Uy{v,z)]dv <G. (A.l) 

Jy 

Since (by Assumption 3.6 ^b)) Uy{v,w) is strictly increasing in the second argument it follows 
that z > z. 



Proof of Proposition 3. 11| Suppose / is u-convex and n-subdifferentiable. Let y be a point of 



differentiability of / and z G d^f{y). Then for /i > 0, we have both f{y+h)—f{y) = f {y)h-\-o{h) 
and 

f{y)- f{y)>u{y,z) -u{y,z) \/y G Dy. 

From the definition of u-convexity /(y) = n(y, z) — f^{z) and f{y + h) > u{y + h, z) — f^{z), so 
taking y = y + h and putting all this together, 

f{y + h) - f{y) > u{y + h,z) - u{y, z) = Uy{y, z) + o{h). 

Dividing by /i > on both sides and letting /i — t- we get 

f'{y) > Uy{y,z). 
If instead we take h < above then we get the reverse inequality and hence 

f'{y) = uy{y,z). (A.2) 

By the Spence-Mirrlees condition Uy{y, z) is injective with respect to z. Hence, the last equation 
determines z uniquely; thus whenever / is differentiable the sub-differential d^f{y) is a singleton 
set. We set 5^/(2/) = {z*{y)}. Then, since z*(y) G ^"/(y) we have /(y)+/"(z*(y)) = u{y,z*{y)) 
as required. 

To prove the converse statement, suppose that for any point of differentiability y of /, 

fiy) = Uyiy,z*iy)) 
where z* non-decreasing. 

Define 

C(y) = sup[/(u) + u{y, z*{v)) - u{v, z*{v))]. 

V 

Then C is u-convex: to see this define 

S = {{z, a); 3w £ Dy \ z = z*{w), a = f{w) — u{w, z*{'w))} 
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and then 

C(y) = sup [u{y,e) + a]. 

{e,a)£S 

Clearly (,> f . We want to show that Q = f ■ By the assumptions on / we have 

ry 
f{y)= Uy{w,z*{w))dw + f{v). 

J V 

Hence 

ry 
f{y) - {f{v) + u{y, z*{v)) - u{v, z*{v))) = / [uy{w, z*{w)) - Uy{w, z*{w))]dw > 0. 

Jv 

Thus C = / &iid so / is n-convex. 



The following Corollary is immediate from 3.11 



Corollary A. 2. Suppose f is differentiable at y. Then f is twice- differentiable at y if and only 
if z* is differentiable at y. 



Proof of Proposition 3.13 



Suppose / is n-subdifferentiable in a neighbourhood of y. Then for small enough e, 

fiy + e) - fiy) > f{y + e, z*{y)) - fiy, z*{y)) 

and lim4o{/(y + e) - f{y)}/e > fy{y,z*{y)). 

For the reverse inequality, if z* is continuous at y then for e small enough that z*{y+e) < z*{y)+5 
we have 

f{y + e)- f{y) < f{y + e, z*{y + e)) - /(y, z*{y + e)) < f{y + e, z*{y) + 5)- f{y, z*{y) + ,5) 

and lim4o{/(y + e) - /(y)}/e < lims^o fy{y,z*{y) + 6) = fy{y,z*{y)). 

Inequalities for the left-derivative follow similarly, and then f'{y) = Uy{y,z*{y)) which is con- 
tinuous. 

Conversely, if d^f is multi- valued at y so that z* is discontinuous at y, then 

lun{f{y + e) - f{y)}/e > fy{y, z*{y)+) > fy{y, z*{y)-) > lim{/(y) - f{y - e)}/e 



where the strict middle inequality follows immediately from Assumption 3.6 
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